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Abstract 

Supersymmetric quantum mechanics (SUSY QM) is a powerful tool for generating 
new potentials with known spectra departing from an initial solvable one. In these 
lecture notes we will present some general formulae concerning SUSY QM of first and 
second order for one-dimensional arbitrary systems, and we will illustrate the method 
through the trigonometric Poschl- Teller potentials. Some intrinsically related subjects, 
as the algebraic structure inherited by the new Hamiltonians and the corresponding 
coherent states will be analyzed. The technique will be as well implemented for peri- 
odic potentials, for which the corresponding spectrum is composed of allowed bands 
separated by energy gaps. 

1 Introduction 

The techniques based on the factorization method, which try to identify the class of Hamil- 
tonians whose spectral problem can be algebraically solved, have attracted for years people's 
attention. In these lecture notes we will elaborate a brief survey of this subject; since we 
shall show that the factorization method, supersymmetric quantum mechanics (SUSY QM), 
and intertwining technique are equivalent procedures for generating new solvable poten- 
tials departing from an initial one, these names will be indistinctly used to characterize the 
technique we are interested in. For a collection of books and review articles concerning 
factorization method and related subjects the reader can seek [1]- [24]. 

I would like to emphasize that these lecture notes on SUSY QM are a sequel of those 
delivered by myself in 2004 at the Latin- American School of Physics [23]. For people wanting 
to have a global view of the way our group addresses the subject, and to know more details of 
our achievements, I recommend to read both, this article and [23], because they complement 
to each other. 

From a historical viewpoint, it is nowadays accepted that the factorization method was 
introduced for the first time by Dirac in 1935 in order to derive algebraically the spectrum 
of the harmonic oscillator [25]. In coordinates such that h = u = m = 1 the oscillator 
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Hamiltonian admits the following factorizations: 
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From these equations it is straightforward to derive the intertwining relationships, 

Ea^a{E-l), Ha+ ^a+{H +1), 
which lead to the commutators ruhng the standard Heisenberg-Weyl algebra: 

[H,a] = —a, [H,a~^] = a"*", [0,0"*"] = 1. 
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These expressions and the fact that H is positive definite allow to generate the eigenfunctions 
ijjn{x) and corresponding eigenvalues En = n -\- 1/2, n = 0,1, . . . oi H . 

The second important advance in the subject was done by Schrodinger in 1940, who 
realized that the factorization method can be as well applied to the radial part of the Coulomb 
problem [26, 27] . In this case the explicit factorizations in appropriate units read: 
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As in the previous case, these factorizations imply some intertwining relationships: 

He+id^'^i = al^^Hi, Hi^iai = aiHi, (8) 

from which it can be derived once again the radial eigenfunctions Rne{r) and eigenvalues 
En = -1/(2^2), n = £+!,£ + 2, .. . for the Hamiltonian iJ^ 

In a set of works starting from 1941, Infeld and collaborators put forward the method 
[28,29], culminating their research with a seminal paper in which it was classified a wide set 
of Schrodinger-type Hamiltonians (indeed four families) solvable through the factorization 
method [1]. The main tool for performing this classification were the factorizations 
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which lead straightforwardly to the intertwining relationships 
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With these ingredients, Infeld and Hull were able to derive elegantly the eigenfunctions and 
eigenvalues for the hierarchy of Hamiltonians Hm [!]• Moreover, the very technique allowed 
them to determine the several forms of the potentials Vm{x) admitting the factorization 
treatment. In this way, after [1] the idea that the factorization methods was exhausted 
started to dominate; indeed, at some time it was thought that for a potential V{x) to be 
solvable through factorization, it just has to appear in the Infeld-HuU classification. 

Reveahng himself against this idea, in 1984 Mielnik performed a simple generahzation of 
the factorization method [30] . He proposed to look for the most general first-order differential 
operators 6, factorizing the harmonic oscillator Hamiltonian in the way: 



H = bb^ 
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Mielnik was able to determine the most general form of /3{x). Moreover, it turns out that 
the product b'^b supplies a new Hamiltonian 
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V{x) = P'{x) + 1. 



(14) 
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Hence, the following intertwining relationships are valid: 

Hb+ ^b+(H + 1), Hb^b{H-l), 



(16) 



which interrelate the eigenfunctions of the harmonic oscillator Hamiltonian H with those of 
H and vice versa. 

Mielnik's work was a breakthrough in the development of the factorization method, 
since it opened new ways to explore exactly solvable potentials in quantum mechanics. In 
particular, his generalization was quickly applied to the radial part of the Coulomb problem 
[31]. The important resultant expressions for that system become: 



Hi = A+A^ - — , 
He-iAi 
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2£2' 



H.Aj = AjH,.,. 



(17) 



They allow to generate the eigenfunctions and eigenvalues of i/^-i departing from those of 
Hf. Notice that, in the same year the connections between the factorization method with the 
Darboux transformation and Witten's supersymmetric quantum mechanics were established 
by Andrianov et al [32, 33] and Nieto [34] respectively. 

In 1985 Sukumar [35,36] pushed further Mielnik's factorization by applying it to general 
one-dimensional potentials V{x) and arbitrary factorization energies e in the way 



H = AA+ + e, A+A + e, 
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These relationships imply that 

HA+ = A+H, HA = AH, (22) 

which allow to generate the eigenf unctions and eigenvalues of H from those of H. 

Along the years, several particular cases and applications of the generalized factorization 
of [30,31,35] have been performed successfully [37]- [63]. Up to this point, however, the 
factorization operators, becoming the intertwiners in the equivalent formalism, were differ- 
ential operators of first order. A natural generalization, posed in 1993 by Andrianov and 
collaborators [64, 65] , requires that Eq. (1221) be valid but A, A~^ are substituted by differential 
operators of order greater than one. Through this higher-order case, it is possible to surpass 
the restriction arising from the first-order method that we can only modify the ground state 
energy of the initial Hamiltonian. In 1995 Bagrov and Samsonov proposed the same general- 
ization, and they offered an alternative view to this issue [66] . On the other hand, our group 
arrived to the subject in 1997 [67]- [71], although there were several previous works related 
somehow to this generalization [72]- [76]. It is worth to mention explicitly some important 
contributions that members of our group have made to the factorization method and related 
subjects. 

• Generation of SUSY partner of the harmonic oscillator [30,68,69] 

• Construction of coherent states for SUSY partners of the harmonic oscillator [73-75,77] 
and for general one- dimensional Hamiltonians [78]. 

• Determination of SUSY partners of the radial oscillator potential [76] 

• Characterization of systems ruled by polynomial Heisenberg algebras and connection 
with SUSY partners of the oscillator [77,79-81] 

• Analysis of SUSY partners of the radial Coulomb problem [31,82-84] 

• Implementation of SUSY transformations involving complex factorization 'energies' 
[85]- [87] 

• Study of the confluent algorithm for second-order SUSY transformations and applica- 
tions [84,88,89] 

• Generation of SUSY partners of the Poschl- Teller potentials [90]- [92] 

• Application of SUSY techniques to the Lame and associated Lame potentials [93]- [99] 

It is important to mention that some other groups have elaborated the same subject from 
different viewpoints, e.g., A^-fold supersymmetry by Tanaka and collaborators [100]- [104], 
hidden nonlinear supersymmetries by Plyushchay et al [105]- [109], etc. 
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2 Standard supersymmetric quantum mechanics 



In this section we are going to show that the factorization method, intertwining technique and 
supersymmetric quantum mechanics are equivalent procedures for generating new solvable 
potentials departing from an initial one. The most straightforward way to see this equiva- 
lence is starting from the technique which involves just first-order differential intertwining 
operators. 



2.1 Intertwining technique 

Suppose we have two Schrodinger Hamiltonians 

which are intertwined by first-order differential operators Ai,A^ in the way 



H,At = At Ho, HoA, = A^H^, 



d . , 
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d . , 
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By employing that, at the operator level, it is valid 

iJiiAj KAJiAj KAJiAj \Aj<Xj \Aj<Xj 



(26) 



where f{x) is a multiplicative operator when acting on a generic wavefunction ip{x), the 
calculation of each member of the first intertwining relationship in Eq. (1241) leads to 



Thus: 



1-^^ 1 d^ ai d"^ ,. d a'.' 



V, = Vo-a[, 



By plugging the V\ of Eq. fl29l) in Eq. fl30|) and integrating the result we get: 

a'l + a? = 2[Vo(a;) -e]. 



(27) 
(28) 
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This nonlinear first-order differential equation for a\ is the well-known Riccati equation. In 
terms of a seed Schrodinger solution u^^^ such that ai = u^^^ /u^^^ it turns out that: 



(0)^,^(0) 



(32) 



which is the initial stationary Schrodinger equation. These equations allow to establish the 
link between the intertwining technique and the factorization method, since it is valid: 

Ho = AiA^ + e, Hi= A^Ai + e. (33) 

Suppose now that Hq is a given solvable Hamiltonian, with known eigenfunctions ipn^ 
and eigenvalues En, namely, 

HoiP^n^ = Enip^^\ n = Q,l,... (34) 

In addition, let us take a nodeless mathematical eigenfunction u^^^ associated to e < Eq in 
order to implement the intertwining procedure. The equations fl24|) and the resulting factor- 
izations ([33]) indicate that, if ^ 0, then [ip^n^ = A+t/'I^VVK - e} is an orthonormal 
set of eigenfunctions of Hi with eigenvalues {E^}- This is a basis if there is not a square- 
integrable function ipe^^ which is orthogonal to the full set. Thus, let us look for tpi^^ such 
that 

= (#\#) oc (^e^'\ A+V^i°)) = iAM'\^l^^) aM'^ = 0. (35) 

By solving this first-order differential equation we get: 

^(1) oc e-^o'"i(^)'^^ = 4^. (36) 

Since Hiipi^^ = eipe^^ (compare Eq. (l33!) ). then the spectrum of Hi depends of the normaliz- 
ability of We can identify three different situations. 

(i) Let us choose e = Eq and = which is nodeless in the domain of Vq, ai = 

i^o'^ /■^o"''- Thus, Vi = Vq — o'i is completely determined, and the corresponding Hamiltonian 
Hi has eigenfunctions and eigenvalues given by 

^« = ^^I^, Sp{Hi) = {En, n=l,2,...}. (37) 

Notice that Eq ^ Sp(ifi) since ipi^^ oc l/u'-^^ is not square-integrable. 

(ii) Let now a nodeless seed solution m*-"-* for e < Eq he chosen so that ai = m^"^ /u^^^ is 
non-singular. Note that u'^^^ diverges at both ends of the x-domain =^ ipe^^ oc I/m^^^ tends to 
zero at those ends. Hence, the full set of eigenfunctions and eigenvalues is 

V'Woc-^, ^(^) = -^^, Sp{Hi) = {e. En, n = 0,1,2,...}. (38) 

(iii) Let us employ a solution u^^^ with a node at one of the ends of the x-domain for t < Eq. 
Thus, the transformation induced by ai = u^^^' /u^^^ is non-singular, ipi^^ oc diverges at 
the end where u^^^ tends to zero, and then the eigenfunctions and eigenvalues of Hi become 

= -j^^, Sp{Hi) = {En, n = 0,1,2,...}. (39) 
vEn - e 

This implies that Hi and Hq are isospectral Hamiltonians. 
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2.2 Supersymmetric quantum mechanics 

Let us build up now the so-called supersymmetric quantum mechanics, by realizing the 
Witten supersymmetry algebra with two generators [110] 



[Qi, ifss] = 0, {Qi, Qj} = SijHss, i,j = 1, 2, 



in the way 







o -91±R o -91zR o- 



Hss = {Q,Q^} 



H+ 
H- 



(A+A, 
V A,Ai 
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Hi-t 
H^-e 



Notice that there is a linear relationship between i/gg and = Diag{ifi, -f^o}: 



{HP 



(40) 

(41) 
(42) 

(43) 



There are two-fold degenerated levels of Hss [En — Eq with n = 1, 2, . . . for the case (i) and 
En — e with n = 0, 1, 2, . . . for cases (ii) and (iii)] which have associated the two orthonormal 
eigenvectors 







(0) > 



(1) 







(44) 



The ground state energy of Hss in cases (i) and (ii) is non-degenerated (and equal to 0), with 
the corresponding eigenstate being given respectively by 







^1 



(0) j ' 



(1) 







(45) 



In this case it is said that the supersymmetry is unbroken. On the other hand, in case (iii) 
the ground state energy {Eq — e) is doubly degenerated, and then the supersymmetry is 
spontaneously broken. 



3 Example: trigonometric Poschl- Teller potentials 

Let us apply the previous technique to the trigonometric Poschl- Teller potentials [91,92]: 

(A-l)A 



2sin (x) 2cos2(a;) 



0<x<7r/2, A,z/>1. 



(46) 



First of all, we require the general solution of the Schrodinger equation (1321) with the Vo(x) 
given in P6|) . which is given by 



sm^{x) cos''(x)|A 2F1 
1 + u-X 



2 ^ V2' 2 



e 1+u-X 
2' 2 



f A 1 • 2/ \ 

^■,X + -;sm [x) 



|; ^ - A; sm\x) 



(47) 



where = X + u. The physical eigenf unctions ipn (x) satisfy ipn (0) = ipn (vr/2) = 0. If we 
ask that ijjn\o) = 0, it turns out that B = 0. Moreover, in order to avoid that the divergent 
behavior at a; = 7r/2 of the hypergeometric function of the remaining term dominates over 
the null one induced by cos'^(a;) for ii^ > it has to happen that 

|±yf=-„ ^ E„=(^i^. „ = 0.1.2.... (48) 

The corresponding normalized eigenf unctions are given by 

sin^(x) cos''(x) 2-^i[-n,n + /x; A + i; sin^(x)]. (49) 

It will be important later to know the number of zeros which u^^^ has, according to the 
value taken by e. This information can be obtained by comparing the asymptotic behavior 
of u^'^^ at X ^ and x 7r/2; by picking up -B = 1 and A = —b/a + g in Eq. (H7I) with 

r(A + i)r(.-i) r(|-A)r(.-i) 

it turns out that: 



(0) 



2(^l+2n)n\T{^l+n){\+^)r, 



• If e < Eq^ u^^^ has either or 1 nodes for g > or g < respectively 

• In general, if -Ei-i < e < Ei, u^^^ will have either i or i + 1 nodes for g > or g < 
respectively, i = 2,3, . . . 

Concerning the spectral modifications which can be induced through the standard first- 
order supersymmetric quantum mechanics, it is obtained the following. 

(a) Deleting the ground state of Hq. Let us choose e = Eq and 

= oc sin^(x) cos''(x). (51) 
The first-order SUSY partner of Vq becomes (see Eq.( l29|) ) 



+ A..>1. (52) 

2sm^(a;) 2cos^(x) 

Since ^ oc 1/4°^ diverges at x = 0,n/2 ^ Eq ^ Sp{Hi) = {E^, n = l,2,...}. Note that 
Vi can be obtained from Vq by changing A— s>A + l, u ^ u + 1; this property is known 
nowadays as shape invariance [111], although it was known since Infeld-Hull work [1]. 

(b) Creating a new ground state. Let us take now e < Eq and a nodeless u'^^^ (Eq.( l471) with 
B = 1,A = -b/a + q, q > 0). As u^^^ oo when x 0, 7r/2 ^^^(0) = -^^{71/2) = 0, i.e., 
ipi^^ is an eigenfunction of Hi with eigenvalue e, and therefore Sp(ifi)={e, En, n = 0,1, . . .}. 
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In order to deal with the singularities at x = 0, 7r/2 induced by u'^^^ on Vi, it is expressed 



u 



(0) 



sin (x) cos ''{x)v, where v is a nodeless bounded function in [0,7r/2]. Hence 

(A-2)(A-1) (.-2)(.-l) 

2sm\x) 2cos2(a;) ^ ^ ' 



(53) 



(c) Isospectral potentials. They arise from the previous case (for e < Eq) in the limit when 
M*^*^-* acquires a node at one of the ends of the domain so that ipi^'^ leaves to be an eigenstate of 
Hi. We can take, for example, the solution given in Eq.l HTl) with A = 1, B = ^ u''^\0) = 0. 
We isolate the singularities induced on Vi by expressing u^^^ 
is a nodeless bounded function in [0,7r/2]. Hence 



sin'^(x)cos^ '^(x)v, where v 



A(A + 1) ^ (u~2)iu-l) 
2sin^(a;) 2cos2(x) 



[Inv]", 



A > l,z/ > 2. 



(54) 



Notice that now Hi and Hq are isospectral. 

An illustration of the several first-order SUSY partner potentials Vi{x) of Vq{x) (see 
Eqs. (152115^ ) is given in Figure 1. 



200 



100 




Figure 1: Trigonometric Poschl- Teller potential Vo{x) for X = 3,h> = A (gray curve) and 
its first-order SUSY partners Vi{x) generated through the following spectral modifications: 
(a) deleting the ground state Eq = 24.5 (dashed curve); (b) creating a new ground state at 
e = 19 (continuous curve); (c) isospectral transformation for e = 19 (dotted curve) [92] 



4 Second-order Supersymmetric Quantum Mechanics 

The second-order supersymmetric quantum mechanics can be implemented through the it- 
eration of two first-order transformations. However, by means of this approach it is not 
revealed in full the rich set of spectral manipulations which is offered by the second-order 
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technique. This is avoided if we address the method directly, by assuming that the inter- 
twining operator B2 is of second-order in the derivatives, namely [67,68,83, 112] 

= B^Ho, (55) 
H, = ~-^ + V,{x), 2 = 0,2, (56) 

i^^4(^-.(.)^ + 7(.)). (57) 
A calculation similar as for the first-order case leads to the following set of equations [23] : 

V2 = Vo-v', (58) 
Y-^'-vV2 = 2V^-vVo, (59) 

7V2 - Y = K," - VV^ + 7K). (60) 
By decoupling this system, assuming that Vq is a given solvable potential, we arrive at: 

^ = !^ + t.2Vo + d, (61) 

^ - ^ + J + r/V - 2K)r/2 + dri'' + c = 0, c,deR. (62) 

Notice that V2 and 7 become completely determined once it is found rj; this task is performed 
by assuming the ansatz 

r^' = -rf + 2/3r/ + 2^, (63) 
where the functions j3 and ^ are to be determined. Therefore 

e^c, /5' + /52 = 2[yo-e], e = (rf + 0/2, (64) 

which is once again a Riccati equation for /?. By assuming that (3 = it turns out 

that u^^^ satisfies the stationary Schrodinger equation 

Since ^ = ±a/c, we get in general two different values for e, ei = (c? + -\/c)/2, 62 = (c? — 
a/c)/2. We thus arrive at a natural classification of the second-order SUSY transformations, 
depending on the sign taken by the parameter c [112]. 



4.1 (i) Real case with c > 

For c > we have that ei, 62 G M, ei 7^ e2; denoting by /32 the corresponding solutions of 
the Riccati equation it turns out that (compare Eq. fl63|) ) 

^' = _^2^2/?ir/ + 2(ei-e2), (66) 
r?' = -r?2 + 2/?2r/ + 2(e2-ei). (67) 
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By making the difference of botli equations it is obtained: 

Notice tliat now tlie new potential V2 has no extra singularities with respect to Vq if 
W{uf'\uf^) is nodeless (compare Eq. (!58!) ). 

The spectrum of H2 depends on the normalizability of the two eigenf unctions ipei'2 of H2 
associated to 61^2 which belong as well to the kernel of B2, 

B2^if = 0, H24f = eM'\ J = 1,2. (69) 
Their explicit expressions in terms of uf^ and u'^^ become [23] 



(0) (0) 
W;(2) oc oc fro) 



We have observed several possible spectral modifications induced on the Poschl- Teller po- 
tentials. 



(a) Deleting 2 neighbor levels. By taking as seeds two physical eigenf unctions of Hq associ- 
ated to a pair of neighbor eigenvalues, namely ei = €2 = -Ei-i, u{'' = '4'f\ '^2°'* = ^i-i 

is straightforward to show that 

W{ui, U2) oc sin2^+i(x) cos2'^+^(x) W, (71) 

where 

^ ^ W{2F^[-i, i + + sin^(x)], 2Fi[-i + 1, ^ - 1 + /i; A + |; sin^(3:)]} 

sin(x) cos(x) 

is a nodeless bounded funcion in [0,7r/2]. Hence, 

(A + l)(A + 2) ^ (.+ l)(. + 2) 

2sm^(a;) 2cos2(x) 

(2^ (2^ 

Notice that the two mathematical eigenfunctions ipli , ip\2 of H2 associated to ei = 
Ei, €2 = Ei_i do not vanish at 0, 7r/2, 

hm K)j = oo. (74) 

Thus, -Ei-i, Ei ^ Sp(if2) = {-£^0, . . . -Ei-2, -Ej+i, . . . }, i.e., we have deleted the levels -Ej-i, Ei 
of ifo in order to generate V2{x). 

(b) Creating two new levels. Let us choose now Ei_i < €2 < ei < Ei,i = 0,1,2,..., 
E^i = —00, with the corresponding seeds u^^^ as given in Eq. lHTI) with B12 = 1, ^1,2 = 
— ^i,2/cti,2 + 91,2, (I2 < 0, gi > 0. This choice ensures that uf^ and u^^^ have i + 1 and i 
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alternating nodes, which produces a nodeless Wronskian [23,96, 113]. The divergence of the 
li^'^) solutions for X — > 0, 7r/2 is isolated by expressing 

u^^^ = sin^-^(x) cos^-"(x)vi,2, (75) 

Vi^2 being bounded for x e [0,7r/2], Vi^2(0) 7^ 7^ Vi^2(7r/2). Since the second term of the 
Taylor series expansion of V12 is proportional to sin^(x) =^ tends to zero as sin(x) for 
X — > and as cos(a;) for x 7r/2. Thus, 

W{uf\uf^) = sm^-^\x) cos3-''^(,t) W, (76) 
where W — l^(vi, V2)/[sin(x) cos(a;)] is nodeless bounded in [0,7r/2]. Hence 

Vi=(A^l)(^+«^_(,„M.)". A,.>3. (77) 

2sin''(x) 2cos2(x) 

Since linix-^o.f V'ci,2 = ^ Sp{H2) = {Eq, . . . , Ei_i, 62, ei, Ei, . . .}, i.e., we have created two 
new levels £1,62 between the originally neighbor ones Ei_i,Ei for generating V2{x). 

It is worth to mention some other spectral modifications which can be straightforwardly 
implemented [92]: 

• Isospectral transformations 

• Creating a new level 

• Moving an arbitrary physical level 

• Deleting an original physical level 

An illustration of the potentials V2{x) induced by some of the previously discussed spec- 
tral modifications is given in Figure 2. 

4.2 (ii) Confluent case with c = 

For ^ = we have that e = ei = £2 G M. It this case, once the solution /3 to the Riccati 
equation for e is found, one has to solve the Bernoulli equation resultant for 77, 

r]'^-r,' + 2Pr, 77 = ^, (78) 

w{x) 



(x) ^Wo + j exp[2 J P{x)dx]dx = ^0 + ^ [w^^Hz/)]" (79) 

[Jp{x)dx], (80) 



u^^^ oc exp[ 
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Figure 2: Trigonometric Poschl- Teller potential Vo{x) for A = 5, z/ = 8 (gray curve) and some 
second-order SUSY partners V2{x) in the real case generated by producing the following 
spectral modifications: (a) deleting the levels E2 = 144.5, E3 = 180.5 (dashed curve); (b) 
creating two new levels at ei = 128, €2 = 115.52 (continuous curve); (e) moving the level 
E2 = 144.5 up to ei = 169.28 (dotted curve) [92] 



with wq G M. Notice that now w{x) has to be nodeless, which is obtained for seed solutions 
u^^^ such that 

limM(°) = 0, / [u^'^\y)]'^ dy < 00 or (81) 

^^^^ Jxo 

lim = 0, / \u^''\y)]^dy < 00, (82) 

Jxi 

and adjusting as well Wq, where [ the x-domain of the involved problem, xq G [xi, Xr]. 

There is one eigenfunction, ipi^^ oc u*^*^) / w, of H2 associated to e belonging as well to the kernel 
of B2, which implies that Sp(if2) depends of the normalizability of ipe ■ 

For the Poschl- Teller potentials, several possibilities for modifying the initial spectrum 
have been observed. 

(a) Creating a new level. Let us choose 3 e ^ and u^^^ as given in Eq. (H71) with 
A = l, B = 0: 

,(0) _ ,l^A/ x^^,^/ N /^/f , ^ /£ _ _ /£. X , 1. .^2^ 



u^'^ = sin\x)cos^{x)2F,i^^ + ^-,^-^-;X + -;sm\x)j (83) 

= sin^(x)cos^"''(x)v(x), (84) 

where v is bounded for x G [0, 7r/2], v(0) 7^ 7^ v(7r/2). By calculating the involved integral 
with xo = it is obtained: 



W = Wq+ ^ ■ 

m=0 



°° ^^+^/f)™(^-^/f)„.sin2^+'"+Hx) 



2^™V 2 

(A + i)„m!(2A + 2m + l) 



^ / l + A — z/ /e l + X — u /e , 11, 3 9, ^ \ , s 

X3F2 ( J-, +W-,A+m+-;A + -,A+m+-;sm2(x) ) . (85) 
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Notice that w is nodeless in [0, 7r/2] for > and it has one node for tuo < 0. By choosing 
now the w of Eq. flHSl) with > 0, we have to isolate anyways a divergence of kind cos 



X 



arising for x 7r/2 through the factorization 

= COs3-2^(x) W, (86) 

W being a nodeless bounded function for x G [0, 7r/2]. From Eg. (1581) we finally get 

(A-l)A (z/-3)(z/-2) 



Vo 



(InW)", A>l,z/>3. 



7) 



2sin^(a;) 2cos'^{x) 

('2') 

Since lim^^^o,- W = 0, it turns out that Sp(if2) = {e, -E-n, n = 0, 1, . . . }. It is important to 



remark that the procedure can be implemented for solutions associated to e > Eq, i.e., by 
means of the confluent second-order transformation one is able to modify the levels above 
the ground state energy of Hq 

Some additional possibilities of spectral manipulation are now just pointed out [92] . 

• Isospectral transformations 

• Deletion of an arbitrary level 

The potentials V2{x) induced by some of the mentioned spectral modifications are illus- 
trated in Figure 3. 




Figure 3: Trigonometric Poschl- Teller potential Vq{x) for A = 5, = 8 (gray curve) and 
its second-order SUSY partners V2(x) in the confluent case, generated through the follow- 
ing spectral modifications: (a) creating a new level at e = 147.92 (continuous curve); (b) 
isospectral transformation with e = 162 (dotted curve); (c) deleting the level = 180.5 
(dashed curve) [92] 
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4.3 (iii) Complex case with c < 

For c < it happens that ei^2 G C, 62 = ci. As we want that V2 to be real, it has to be 
chosen P2 = Pi. Thus, 

2Im(ei)_«;' ^_W{uf\uf^) 



Im[/5i] w' 2(ei-ei) ' ^ ' 

where uf'^ is a complex solution of the Schrodinger equation associated to ei G C As in the 
previous cases, w has to be nodeless in order to avoid the arising of extra singularities in V2. 
For this to happens, it is sufficient that 

lim uf^ = or lim wf^ = 0. (89) 

In both cases we get that Sp{H2) = Sp(ifo)- 

We illustrate now the technique through the trigonometric Poschl- Teller potentials, for 
which we take e G C, as in Eq.(^ with A = l, B = ^ mW(0) = 0, 00 when 

X ^ n/2. It is appropriate to express now 

= sin^(x) cos^-''(x)v. (90) 

Hence, it is obtained: 

w = sin2^+^(x) cos3-2'^(a;) W, (91) 

W = ^^^'^^ , (92) 

2(e — e) sin(x) cos(x) 

>,.(™±^+(^i_i?l^-„„>V)", A>1..>3. (93) 

2sm (x) 2cos^(x) 

The potentials Vq and V2 turn out to be isospectral. An illustration of Vo(x) and the V2(x) 
of Eq.( l93i) is given in Figure 4. 

A general conclusion, which can be inferred of the previously discussed first and second- 
order SUSY QM but it is valid for an arbitrary order, is that the spectra of the initial 
and new Hamiltonians differ little. Thus, one would expect that some properties which are 
somehow related to the corresponding spectral problems would be essentially the same. We 
will see next that this is true, in particular, for the algebraic structure ruling the two involved 
Hamiltonians. 

5 Algebra of the SUSY partner Hamiltonians 

Given an initial Hamiltonian Hq, which is characterized by a certain algebraic structure, it 
is natural to pose the following questions: 

• What is the algebra ruling the potentials generated through SUSY QM? 

• What other properties of Hq are inherited by its SUSY partner Hamiltonians? 
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-500 L 

Figure 4: Trigonometric Poschl- Teller potential Vo{x) for A = 5, z/ = 8 (gray curve) and its 
second-order SUSY partner V2{x) in the complex case generated by using e = 176.344 + 1.5i 
and a solution m^^^ vanishing at the origin (dotted curve) [92] 



These questions will be analyzed here through a simple system which algebraic structure 
generalizes the standard Heisenberg-Weyl (oscillator) algebra. Thus, let us consider a one- 
dimensional Schrodinger Hamiltonian 

Ho = -l^,+Vo{x), (94) 

whose eigenvectors {ipn^) and eigenvalues En satisfy (for calculational convenience we switch 
now to the Dirac bra and ket notation) 

i/o|V'i°^) =i?n|^i°^ Eo<E,<E2<... (95) 

oo 

= Smn, Yl = 1- (96) 

m=0 

It is supposed that there is an analytic dependence between the eigenvalues En and the index 
labeling them, namely. 

En = E{n). (97) 

Let us define now the intrinsic algebra of Hq, by introducing a pair of annihilation and 
creation operators whose action onto the eigenvector of Hq is given by 

ao l^f ) = r7(n)|^i°i), ao+| ^l°))=f,(n + l)|^i°i), (98) 
rj{n) = e^°[^W-^("-i)l ^E{n) - Eq, a G M. (99) 

With these definitions it is straightforward to show that: 

a+ao=Ho-Eo. (100) 
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It is introduced the number operator Nq through its action onto the eigenvectors of Hq: 

TVolV^f ) = nlV^f ), [No,4] = ±4- (101) 

Now, the intrinsic algebra of the systems is characterized by 

a+oo = E{No) - Eo, a+ = E{No + 1) - Eo, (102) 
K- , a+] = E(No + 1) - E(No) = f{No). (103) 

It is important to realize that for general E{n) the intrinsic algebra becomes nonlinear (since 
f{No) is not necessarily a linear function of A^o)- However, this intrinsic nonlinear algebra 
can be linearized by deforming the corresponding annihilation and creation operators in the 
way 

_ rLjNo + 1) _ r^iVo) _ 

~ rj{No + 1) ° ' ~ rj(iVo) ° ' ~ ^ ' ^ ^ 

[No, 4j = ±a^^, a+ = No, = Nq + 1, [a^^, a+ ] = 1. (105) 

In particular, for the harmonic oscillator Hamiltonian the intrinsic algebraic structure 
reduces to the standard Heisenberg-Weyl algebra: 

E{n)^n+l-, f{n) ^ E{n + 1) - E{n) ^1, (106) 

[No,a^] = ±a^, [ao,a+] = l. (107) 

On the other hand, for the trigonometric Poschl- Teller potentials we recover the su(l,l) 
algebra: 

E{n) ^ (^+^^)' , f{n)^E{n + l)-E{n)^4:n + 2fx + 2, (108) 

[TVo, 4] = ±4, [«o", 4] = 47Vo + 2// + 2. (109) 

Let us characterize now the algebraic structure of the SUSY partner Hamiltonians of 
ifo [78]. In order to illustrate the procedure, we restrict ourselves just to a second-order 
SUSY transformation which creates two new levels at the positions ei, 62 for H^, namely, 

Sp(//2) = {ei, 62, n = 0, 1, . . . } (110) 

Notice that the eigenstates of H2, l^/'n ), associated to En,n = 0, 1, . . . are obtained from 
[■^i^^) and vice versa through the action of the intertwining operators , B2 in the way: 

WW) = u(o,) = ^ (111) 

The eigenstates of H2 constitute a complete orthonormal set, 

(V'?^lV'i'^) = 0, = {i^fMf^)^5,„ n,m = 0,l,..., i,j = l,2, (113) 

2 oo 
1=1 m=0 
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It is important to notice that the two new levels ei, €2 can be placed at positions essentially 
arbitrary (either both below the initial ground state Eq or in between two neighbor physical 
levels Ei_i,Ei). This fact somehow breaks the symmetry defined by E{n). Thus, it should 
be clear that we need to isolate, in a sense, the two new levels €1,62. With this aim, it is 
defined now the number operator N2 through the following action onto the eigenstates of 
H2: 



iV2|V^f)=n|^(2)), iV2|V^(f)=0, n 



0,1,... 2 



1,2. 



(115) 



Departing now from the operators of the intrinsic algebra of Hq, let us construct those of 
the natural algebra of H2 [78], which represent a generalization of the operators introduced 
previously for the SUSY partners of the harmonic oscillator [30,73,74,77], 



Btar. B2. 



(116) 



An schematic representation of these operators is given in Figure 5. Their action onto the 
eigenvectors of H2 is given by 



Ho 



at 



{2): 

I {2)'. 

m. ) 






< 


— ► 














< 


— > 





a 







(0)^ 



Figure 5: Representation of the annihilation and creation operators of the natural algebra 
of Ho. 



= 0, 02,l^f ) = rN{n) a+^l^i^)) = f^{n + 1) (117) 



rN[n) 



\{[E{n) - ti][E{n - I) - ti] \ r,(n), z = l,2, n = 0,l. 



:il8) 



.1=1 
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Using these equations it is found that 

oo 

= [r>(iV2+l)r^(iV2+l) -f^(iV2)r^.(iV2)] J]|^^2)X^^2) |. (119) 



m=0 



Let us build up now the annihilation and creation operators af which generate the 
intrinsic algebra of H21 



2) 



rN{N2 + 1) rjv(A^2 



[a2,at]^f{N2)J2\^^r^^)ii^^^\ (120) 



m=0 



Their action onto the eigenstates of H2 reads 

atl^S^) = 0, l^i^)) = r,(n)|^2,), (121) 

^=l,2, n = 0,l,... (122) 

Notice that this algebra of H2 is the same as the intrinsic algebra of Hq when we work on the 
restriction to the subspace spanned by the eigenvectors iV'n ^) associated to the isospectral 
part of the spectrum. 

Let us construct, finally, the annihilation and creation operators af^ of the linear algebra 
oiH2, 

_ tl{N2 + 1) _ + r^m + 

ri{N2 + l) ^ ri{N2) ^ ^ ' 

a£] = ±a£ , k" , <] = E 1^^'^) (^^'^ I ■ (124) 

m=0 

Their action onto the eigenstates of H2 is given by 

= 0, «2,I^F^) = rL(n)|^f2i), (125) 
4jV'f)=^L(^+ 1)1^11), ^ = l,2, n = 0,l,... (126) 

Once again, this algebra of H2 coincides with the corresponding one of Hq on the subspace 
spanned by {j^'i^^}, n = 0, 1, . . . }. 



6 Coherent states 

Once identified the algebraic structures associated to Hq and i/2, it would be important to 
generate the corresponding coherent states. It is well know that there are several definitions 
of coherent states [114-117]; here we will derive them as eigenstates with complex eigenvalues 
of the annihilation operators of the systems. Since there are available several annihilation op- 
erators for Hq and H21 we will call them with the same name as the corresponding associated 
algebra. 
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6.1 Coherent states of Hq 

Since for Hq we identified two different algebraic structures, intrinsic and linear one, we will 
look for the two corresponding families of coherent states. 

6.1.1 Intrinsic coherent states \z,a)Q of Hq 

The intrinsic coherent states \z,a)Q of Hq satisfy: 

aQ\z,a)Q = z\z,a)Q, zeC (127) 
Expanding \z, a;)o in the basis of eigenstates of Hq, 



oo 
n=0 



)o = X]c.|^W), (128) 



and substituting in Eq. 0127p . it is obtained a recurrence relationship for Cn, 

^ ^-ia{E„-Eo) 

rj{n) ri{n)...rj{l) - ^o) • • • (^i - ^o) 

By choosing cq G so that \z, a)o is normalized we finally get 



\z, a 



)o = I E ^ E e-^^-'^^^m, (130) 



n P"l In V Pi'' 

<^m=Q / m=0 * 



_ / 1 for m = 0, , . 

^"^ \(E„-Eo)...(i5;i-i?o) form>0. ^'^'^ 

An important property that should be satisfied is the so-called completeness relationship, 
namely, 



\z,a)o Q{z,a\dn{z) = 1, (132) 
where the positive defined measure d^{z) can be expressed as 

/ oo : : \ 

d^{z) = -[^^-^\p{\z\^)d''z. (133) 

Working in polar coordinates and changing variables y = l^p, it turns out that p{y) must 
satisfy 

/"OO 

/ y"'piy)dy = p^, m = 0, 1, . . . (134) 
Jo 

This is a moment problem depending on p„, i.e., on the spectral function E{n) [73,77,118- 
122]. If one is able to solve Eg. (11341) for p{y), the completeness relation becomes valid. 
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Once Eg. (11321) is satisfied, an arbitrary quantum state can be expressed in terms of 
coherent states. In particular, an intrinsic coherent state |z',a)o admits the decomposition, 

\z',a)o= / \z,a)o o{z,a\z',a)odfi{z), (135) 



where the reproducing kernel o{z, a\z', a)o is given by 

Vm=0 ^™ / \m=0 ^™ / \m=0 ^™ / 

As \z = 0,a)o = li'o''^), it turns out that the eigenvalue z = is non-degenerate. If the 
system is initially in an intrinsic coherent state \z,a)o, then it will be at any time in an 
intrinsic coherent state, 

f/o(t)|^,«)o = e-**^"|^,« + t)o, (137) 
where Uo{t) = exp{—iHot) is the evolution operator of the system. 

6.1.2 Linear coherent states 1^;, 0)0^ of Hq 

Now, let us look for the linear coherent states such that 

aQ^\z,a)oi^ = z\z,a)o^, zeC. (138) 
A procedure similar to the previous one leads straightforwardly to 

k,«)o. = e-^f;e— (^™-^°)^|^^°)), (139) 

which, up to the phases, has the form of the standard coherent states. Now the completeness 
relationship is automatically satisfied, 

1 



TC 



z,a)Q^ o^{z,a\d z = 1. (140) 



Therefore, an arbitrary linear coherent state \z', q;)o^ is expressed in terms of linear coherent 
states, 

\z\ a)o^ ^\ j "-^^^ °^ "-^"^ ^^"^^^ 

where the reproducing kernel is given by 

Q^{z,a\z' ,q)o^ = exp (-|z|V2 + zz' - \z'\'^/2) . (142) 

Notice that, once again, the only eigenstate of Hq which is also a linear coherent state is the 
ground state, \z = 0,q;)o^ = {ipo^^)- Moreover, since [aot,aQ^] = 1 then the linear coherent 
states also result from the action of the displacement operator Dl{z) onto the ground state 

l4°'>: 

k,a)oi = Dl{z)\iP^°^) = exp(2;a+ - zaQjlip^^^). (143) 
Finally, an initial linear coherent state evolves in time as a linear coherent state, namely, 

Uo{t)\z,a)o, = e-''^'\z,a + t)Q,. (144) 
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6.2 Coherent states of H2 



Let us remember that we have analyzed three different algebraic structures for H2, namely, 
natural, intrinsic and linear ones; thus three families of coherent states will be next con- 
structed. 



6.2.1 Natural coherent states |z,a)2^ of H2 

The natural coherent states \z,a)2^ of H2 are defined by 



A procedure similar to the previously used leads to 



2jv ■ 



(145) 



2jv 



m=0 



1 2m 



m=0 



-ia{Em-E(i) _ 



(146) 



where 



for m = 0, 



Pm fl(^m - ti){Em^i - tif . . .(^1 - 6^)2(^0 " £*) for m > 0. 



(147) 



Notice that now the completeness relationship has to be modified to include the projector 
onto the subspace spanned by the eigenvectors of H2 associated to the two isolated levels 
ei, e2, namely. 



2 r 

/ k,«)2^2^(2;,«|rfM^) = l, 

j=l 



:i48) 



where 



djl{z) = - V -z— 



dy = Pm, m> 0. 



(149) 
(150) 



It is clear that the moment problem characterized by Eg. (11 50 P is more complicated than 
the one associated to Eg. (11341) (compare the moments given by Egs. (11311) and (11471) ). Since 

(2) (2} (2) 

B2\ipei ) = a2]vl'^fi ) = i = 1,2 and a2]vlV'o ) = O5 then the eigenvalue z = of a^^ is 
3-fold degenerate. Finally, as for the coherent states of Hq, it turns out that the natural 
coherent states \z,a)2f^ evolve coherently with t, 



U2{t)\z,a)2, 



-itEo 



z,a + t)2, 



:i5ii 



where U2{t) = exp{—iH2t) is the evolution operator of the system characterized by H2. 
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6.2.2 Intrinsic coherent states \z,a)2 of H2 

The intrinsic coherent states \z^a)2 of H2 are eigcnstatcs with complex eigenvalue 2; G C of 
the intrinsic annihilation operator oi H2. A treatment similar to the previous one leads 
to 

/ °° U|2m\ "2 00 rn 

\ „ rm In V Pm 

\m=0 / m=0 * 

The modified completeness relationship reads now 



+ / \Z:C.)22{z,a\di^{z) = 1, (153) 

1=1 

(2) ('2) 

where d//(2;) is the same as for Hq. Since a2 l^e/ ) =0, i = 1, 2 and \z — 0, q;)2 = \ipQ ), the 
eigenvalue 2; = of 03 turns out to be 3-fold degenerate. As for the natural coherent states 
of i?2, the intrinsic coherent states \z,a)2 evolve in time in a coherent way. 

6.2.3 Linear coherent states \z,a)2i^ of H2 

Let us consider, finally, the linear coherent states \z^a)2^ of H2 as eigenstates of aj^ with 
eigenvalue 2; e C. A standard procedure produces now 

\z, a)2, = e-^ e--(^-^°)4= m. (154) 

The completeness relationship is expressed as 

Y.\^f^){^f^ + lj \z,a)2,2A^,a\d'z = l. (155) 

Once again, the eigenvalue 2; = of a^^ is 3-fold degenerate. Moreover, the linear coherent 
states \z, Ci)2i^ appear from the action of the displacement operator D2^{z) onto the extremal 
state I ■00^'*): 

\z, a)2^ = D2^^{z)\^|J^o^) = exp(^a^^ - 5a2i)lV'S^^)- (156) 
As in the previous cases, the linear coherent states of H2 evolve in time coherently. 

Let us end up these two sections, about the algebraic structure of the SUSY partner 
Hamiltonians and the corresponding coherent states, with the following conclusions. 

• The intrinsic and linear algebraic structures of Hq are inherited by H2 on the subspace 
associated to the levels of the isospectral part of the spectrum 

• This refiects as well in the corresponding coherent states of H2, linear and intrinsic 
ones, whose expressions on the same subspace are equal to the corresponding ones of 
Ho 

• We generahzed successfully the construction of the natural algebra of H2 when the 
initial potential Vo(a:) is more general than the harmonic oscillator 
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7 SUSY partners of periodic potentials 



The exactly solvable models are important since the relevant physical information is encoded 
in a few expressions. Moreover, they are ideal either to test the accuracy (and convergence) of 
numerical techniques, or to implement some analytic approximate methods, as perturbation 
theory. Unfortunately, the number of exactly solvable periodic potentials is small, so it is 
important to enlarge this class. As we saw in the previous sections, a simple method to do 
the job is supersymmetric quantum mechanics. Before applying the technique to periodic 
potentials, however, let us review first the way in which it is classified the energy axis in 
allowed and forbidden energy bands. 

For one-dimensional periodic potentials such that Vo{x + T) — Vo{x), the Schrodinger 
equation 

~r{x) + Vo{x)ilj{x) = E^{x), (157) 



^^>{x) = A{x)^{xl ^>{x) = ( ) , (158) 



can be conveniently expressed in matrix form 

^^^>{x)=A{x)nx,. -W=^^,(^ 

1 

2[Vo{x)-E] 

Therefore: 

*(x) =6(x)*(0), (160) 

where the 2x2 transfer matrix h{x) is symplectic, which implies that its determinant is 
equal to 1. The transfer matrix can be expressed in terms of two real solutions Yi^2{x) such 
that vi(0) = 1, vi(0) = 0, V2(0) = 0, v^(0) = 1, as follows 

- ( v:S ) ^ 

The general behavior of and the kind of spectrum of the Hamiltonian H^^ depends on 
the eigenvalues /3± of the Floquet matrix b{T), which in turn are determined by the so-called 
discriminant D — D{E) — Tr[6(T)]: 



(3^-D(5 + l^Q P^^^±^^-l, /3+/3_ = l. (162) 

By picking up ^E'(O) as one of the eigenvectors of the Floquet matrix b{T) with eigenvalue /3, 
it turns out that 

*(T) =6(T)^(0) =/5^(0). (163) 

Moreover, for these vectors in general we have that 

^{x + nT) = /3''^{x). (164) 

The corresponding t/;(,T) are known as the Bloch functions, and they exist for any G ffi. 
Depending on the value of the discriminant, the Bloch functions acquire three different 
asymptotic behaviors: 
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• If \D{E)\ < 2 for a given E, it turns out that 

(3± = exp(±ifcT), (165) 

where /c G M is called the quasimomentum. Equation (11641) implies that the Bloch 
solutions ip'^i^x) in this case are bounded for \x\ oo, and hence they can be equipped 
with a physical meaning. Therefore, E belong to an allowed energy band. 

• The equation \D{E)\ = 2 defines the band edges, denoted 

Eq < Ei< Ey < E2< E2> < . . . < Ej <Ey<... (166) 

For these energy values it turns out that /3+ = /5_ = ±1, the two Bloch functions tend 
just to one which becomes periodic or antiperiodic. Thus, these energy values belong 
as well to the spectrum of Hq. Moreover, the Bloch eigenfunctions ipj, ipf associated 
to Ej, Eji are real, they have the same number j of nodes, and both are either periodic 
or antiperiodic 

• If \D{E)\ > 2 for a given E, it turns out that 

P+ = P^OeR, p. = ^. (167) 

Then, the Bloch solutions ip^ are unbounded either for x — > oo or for x — — oo so 
that they cannot be equipped with a physical interpretation. Therefore, E belongs to 
a forbidden energy gap 

D(E) 
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Figure 6: The discriminant D{E) versus E for Vq{x) = 5sin^(x) 

As an illustration, we have plotted in Figure 6 the discriminant D{E) as a function of E 
for Vo{x) = 5sin^(x). 
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7.1 Lame potentials 

A very special class of periodic systems are given by the Lame potentials 

Vq{x) = -n{n + l)msii'^{x\m), n&N, (168) 

where sn(a:|m) is a Jacobi elliptic function. These are doubly periodic special functions 
defined by [99] 

sn(a;|m) = sin(^, cn(a;|m) = cos(^, dn(a;|m) = (169) 

ax 

with X = / . The real and imaginary periods are respectively {AK, 2iK'), {AK, AiK'), 

y l—mshP 6 

{2K,4:iK'), where 

7/1 

K = K{m) = I ' =, K' = K{mi), rm^l-m. (170) 

Jo V 1 — m sin^ 9 

The special nature of the Lame potentials is due to the following properties: 

• They have 2n + 1 band edges, which define n + 1 allowed energy bands (n of them 
finite and 1 infinite) and n + 1 energy gaps (n finite and 1 infinite) 

• It is possible to determine the general solution of the Schrodinger equation for any real 
value of the energy parameter 

• In particular, there exist explicit analytic expressions for the band edge eigenf unctions 
Let us illustrate this behavior through the simplest example. 

7.2 The Lame potential with n = 1 

In this case the Bloch functions for a real arbitrary energy value e are given by 

^.(°)^(a;)oc^^^L±i^e-<(«), ^(°)- (x) oc + 7^ e<(°), (171) 

^ _ g2aC(w)-2a;C(a)^ (^172) 

where C^'iz) = —p{z), [In a{z)]' = Ci^): oj — K, u' — iK' and 

2e = ^(m + l)-p(a). (173) 

Note that p is the Weierstrass p-function, which is related to the Jacobi elliptic function sn 
through [123] 

p{z)^es+ ''^^y (174) 
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There are now three band edges, and the exphcit expression for the corresponding eigenvalues 
and eigenfunctions are given by 

Eo = j, ^r(^) = dna;, (175) 

Ei = l, ^f(x) = cna:, (176) 

= i^, ijf!\x) = snx. (177) 

The energy axis is classified as follows: the energy interval (— oo, m/2) represents the infinite 
energy gap (forbidden energies), then it comes the finite energy band [m/2, 1/2] (allowed 
energies), then it arises the finite energy gap (1/2, 1/2 + m/2), and finally it emerges the 
infinite energy band [1/2 + m/2, oo). 



7.2.1 First-order SUSY using Bloch functions 

Let us choose in the first place as seed one real nodeless Bloch function u^^^ with e < Eq. It is 
straightforward to check that u^^^ /u^^^ is periodic, which implies that Vi{x) is also periodic 
(see Eq.( !29|) ). The SUSY transformation maps bounded (unbounded) eigenfunctions of 
Hq into bounded (unbounded) ones of Hi. Since for e < Eq the chosen Bloch function 
y(o) diverges either for x ^ oo or x — oo (e ^ Sp{Hq)), then will diverge at 

the opposite limits and therefore e Sp(iJi). On the other hand, for e = Eq E Sp{Ho), 
the corresponding Bloch solution u^'^^ is nodeless periodic and therefore is also a 

nodeless periodic eigenfunction of Hi with eigenvalue Eq. In conclusion, for first-order SUSY 
transformations involving a nodeless Bloch function with e < E'o it turns out that Sp(ifi) = 
Sp{Hq), i.e., the SUSY mapping is isospectral. 

In particular, for the Lame potential with n = 1 the first-order SUSY transformation 
which employs the lowest band edge eigenfunction ?/'q°''(x) = dnx leads to [124,125] 

Viix) = Vo(^x + ^y (178) 

and it is said that Vi{x) is self-isospectral to Vo{x). This property was discovered, for the first 
time, by Braden and McFarlane [124], although the name is due to Dunne and Feinberg [125]. 
Moreover, by using as seed the Bloch eigenfunction with e < Eo'is is straightforward 

to show that [95] 

Vi{x) = Vo{x±a), (179) 
where a is related to e through Eq. (11731) . 



7.2.2 First-order SUSY using general solutions 

Let us choose now a real nodeless linear combination u^^^ of the two Bloch functions for 
e < Eq. Although now u^^^ /u^^^ is not strictly periodic, it is asymptotically periodic when 
|x| — > oo. This property is transferred as well to the new potential Vi{x) which is asymptoti- 
cally periodic for \x\ oo. As in the previous case, the bounded (unbounded) eigenfunctions 
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of Ho are mapped into bounded (unbounded) eigenfunctions of Hi. Since lim|a;|^oo \u^^^\ = oo, 
it turns out that is square-integrable in (— cxd,cxd). We conclude that, when appro- 

priate hnear combinations of two Bloch functions of Hq associated to e < Eq are used 
the implement the first-order SUSY transformation, we create a bound state for Hi, i.e., 

Sp(ifi)={e}USp(i7o). 

We apply now this kind of transformation to the Lame potential with n = 1, using as seed 
a real nodeless linear combination of the two Bloch solutions u^^^ (x) given in Eq. fll71l) with 
e < Eq [95,96]. In this way it is obtained an asympotically periodic potential, as illustrated 
in Figure 7, where the corresponding bound state of Hi is as well shown. 
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Figure 7: Up. First-order SUSY partner Vi (black curve) of the Lame potential (gray curve) 
for m = 0.5, e = 0.22. Down: The corresponding (unnormalized) bound state of Hi 
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7.2.3 Second-order SUSY using Bloch functions 

Let us apply now a second-order SUSY transformation using two real Bloch solutions uf'l 
associated to ei^2 ^ {—oo,Eq] or ei^2 ^ [Ej^Ej'j^j = 1,2,... so that W{uf'\u^2^) is node- 
less [93]- [96]. For this choice it turns out that W'{uf'\uf^)/W{uf'\uf^) as well as V2{x) 
are periodic (see Eqs. fl58|68p ). The SUSY transformation maps bounded (unbounded) eigen- 
functions of Hq into bounded (unbounded) ones of H2. Thus, Hq and H2 have the same band 
structure, namely, Sp(if2)=Sp(-f^o)- 

For the Lame potential with n = 1, if we employ the two band edge eigenf unctions 
ilj^f\x) , tljf!\x) bounding the finite gap, we get once again [93] 

V2{x) =Vo(^x + ^y (180) 

Moreover, by using as seeds two Bloch functions uf^ (x) such that ei^2 ^ (Ei, Ei') it is 
straightforward to show that [95] 

V2{x) = Vo{x + 6 + 6'), ai = 6 + uj', a2 = 6' + uj'. (181) 



7.2.4 Second-order SUSY using general solutions 

Let us choose now as seeds two real linear combinations of Bloch function for ei,2 G (—00, Eq] 
or ^1,2 £ [Ej,Eji],j = 1,2, . . . so that W{uf'\u':2^) is nodeless [95,96]. Once again, 
W'{uf'\u^^)/W{uf'\uf^) and consequently V2{x) are asymptotically periodic for |x| 00, 
the second-order SUSY transformation maps bounded (unbounded) eigenf unctions of Hq into 
bounded (unbounded) ones of H2. Since now the two eigenf unctions of H2, u^2^ /W{uf\ uf^), 
u^i^ /W{uf'\uf^) associated to ei, €2, can be made square-integrable for x G (—00,00), in 
turns out that Sp{H2)= {ei, 62} U Sp(ifo)- 

Let us illustrate the procedure for the Lame potential with n = 1, using as seeds two real 
linear combinations uf'l of ufl and ufl for €1^2 £ {Ei,Ei'). As in the general case, it is 
obtained a potential V2{x) which is asymptotically periodic for —>■ 00. The corresponding 
Hamiltonian has two bound states at the positions ei, 62. An illustration of these results is 
shown in Figure 8. 

Our conclusions concerning SUSY transformations applied to periodic potentials are the 
following. 

• It was employed successfully the SUSY QM for generating new solvable potentials 
departing from a periodic initial one 

• For seed solutions chosen as Bloch functions associated to the band edges or to factor- 
ization energies inside the spectral gaps, the new potential is periodic, isospectral to 
the initial one 

• For seed solutions which are real linear combinations of Bloch functions, the new po- 
tential is asymptotically periodic, and the spectrum of the corresponding Hamiltonian 
will include a finite number of bound states embedded into the gaps 
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Figure 8: Up. Second-order SUSY partner V2 (black curve) of the Lame potential (gray 
curve) for m — 0.5, ei = 0.64, €2 — 0.65. Down: The (unnormalized) bound state of H2 
associated to €2 — 0.65. 

• The technique was clearly illustrated by means of the Lame potential with n = 1 



8 Conclusions 

• Along this paper it has been shown that supersymmetric quantum mechanics is a 
simple but powerful tool for generating potentials with known spectra departing from 
a given initial solvable one. Indeed, SUSY QM can be used to implement the spectral 
manipulation in quantum mechanics 

• One of our aims when writing this paper is to present the development of the factor- 
ization method, in particular the further advances arising after Mielnik's paper about 
generalized factorizations [30] 
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• Supersymmetric quantum mechanics has interrelations with several interesting areas of 
mathematical physics such as solutions of non-linear differential equations [79-81,126- 
130], deformation of Lie algebras and coherent states [73-75,77,78,131,132], among 
others. In my opinion, a further exploration of these interrelations is needed, which 
represents the future of SUSY QM. We will pursue these issues in the next years with 
the hope of getting as many interesting results as we were obtaining previously. 

Acknowledgments. The author acknowledges Conacyt, project No. 49253-F. 
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